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ABSTRACT 

We consider the RSOS /S-matrices of the *J'{i,5) perturbed minimal models which 
have recently been found in the companion paper These S'-matrices have 
some interesting properties, in particular, unitarity may be broken in a stronger 
sense than seen before, while one of the three classes of *l*(i,5) perturbations (to 
be described) shares the same Thermodynamic Bethe Ansatz as a related ^^(1,2) 
perturbation. We test these new S'-matrices by the standard Truncated Conformal 
Space method, and further observe that in some cases the BA equations for two 
particle energy levels may be continued to complex rapidity to describe (a) single 
particle excitations and (b) complex eigenvalues of the Hamiltonian corresponding 
to non-unitary S-matrix elements. We make some comments on identities between 
characters in the two related models following from the fact that the two perturbed 
theories share the same breather sector. 
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1 Introduction 



This paper is the companion to in which S-matrices were conjectured for ^(1,5) pertur- 

(2) 

bations of Virasoro minimal models by RSOS reduction of the S'-matrix of 03 affine Toda 
theory. Throughout this paper we shall refer to equations in [Q] as (l:nn). 

The description of perturbed conformal field theories (PCFTs) as affine Toda theories 
has provided much insight into their structure — the construction of the S-matrices for the 
particles in PCFTs relies upon the existence of an affine quantum symmetry of the Toda action 

P] — but there are still many interesting aspects to investigate. In this paper we consider 

(2) 

the 02 afhne Toda theory at imaginary coupling which can provide a model for both '^(1,2) 
and ^(1,5) perturbations of minimal models. In the simplest case, the ^(1,2) perturbation of 
the unitary Ising model (which has a minimal eg S'-matrix) is described by the same Toda 
Lagrangian as the $(1,5) perturbation of the non-unitary M^^iq model. 

It is possible for the same Lagrangian to describe two such different models because the 
S-matrices of a PCFT are obtained by 'RSOS restriction' of the Toda S-matrices which relies 
upon & Uq{sl{2)) subalgebra of the quantum symmetry algebra, and that one may be able to 
choose more than one such subalgebra leading to different S-matrices for the PCFTs. In this 
way it turns out that the same Lagrangian theory may indeed describe both a unitary and a 
non-unitary theory, and that even the particle spectra and ground-state thermodynamics of 
these two theories may be the same, yet the theories possess different 5-matrices both being 
consistent solutions of the 5-matrix bootstrap equations for the particle masses. 

(2) 

In the specific case of the 02 Toda theory the RSOS restriction leading to the ^(1,2) PCFTs 
has been performed in ^, ^ and that leading to the $(1,5) perturbation in The way 
that the two perturbations which arise from them are related falls into one of three classes; 
the example considered by Martins in ||^] is in the first class and we give the first examples of 
the other two classes in this paper, together with a unified description. 

(2) 

The 02 affine Toda theory appears to describe a non-unitary theory, and while it is possible 
for the RSOS reduction of a non-unitary theory to yield a unitary scattering theory (as with 
the <&(i,2) perturbations of unitary minimal models) for most reductions the scattering theory 
is also non-unitary. However, for all the reductions proposed so far, this has been a rather 
weak form of non-unitarity, with the Hamiltonian having real eigenvalues bounded below, 
and simply the Hilbert space having an indefinite metric. However, in some of our models, 
we show how unitarity is violated in a stronger sense, and the Hamiltonian of the reduced 
theory may have complex eigenvalues as well. 

The outline of the paper is as follows: in || we describe briefly the minimal models of 
conformal field theory and their perturbations and the way in which two different perturbed 
conformal minimal models, one perturbed by the ^(1,2) field and one by the $(1,5) field, may 
be related to the same affine Toda Lagrangian field theory, and how the BRS reduction of 
the two different underlying Liouville theories naturally leads to two different S-matrices for 
the two perturbed minimal models. 

We then pass to tests of the predicted S-matrices in the case of the ^(1,5) perturbations. In 
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section |^ we recall the predictions of the Thermodynamic Bethe Ansatz for the ground state 
energy of the particle system on a circle, which leads to predictions for the effective central 
charge, the mass-gaps and of the coefficient in the linear term in the asymptotic behaviour of 
the ground-state itself. We also show how in many cases the two perturbed minimal models 
which can be derived from Smirnov's S'-matrix give the same TBA system, while in other 
cases this is not so. We also recall how the higher multi-particle energy levels can be predicted 
using the ordinary Bethe ansatz. 

The first prediction can be checked analytically while the other three can be checked us- 
ing the Truncated Conformal Space Approach (TCSA) of Yurov and Al. Zamolodchikov Q. 
In section ^ we use this method to show that our predictions of the particle spectra and 
S'-matrices are correct in three models, the minimal models ^^3,14 and TWs^ie per- 

turbed by the field *&(i,5)- 

Furthermore we show that the Bethe Ansatz can be usefully continued to complex rapidities 
in two different cases. Firstly, we show that the ordinary Bethe Ansatz may also be a good 
description of a single particle state, when continued to imaginary rapidity. This presumably 
indicates that the Bethe Ansatz equations provide the dominant corrections to the full TBA 
equations for single particle states developed in Q . Secondly, we show that for the non-unitary 
S'-matrices the Bethe Ansatz still provides a good description of the complex eigenvalues 
analogous to two-particle states found in the TCSA investigation. 

In section ^ we give some comments on the possible ways in which the S'-matrix may be 
realised and the connection with the different partition functions and periodicity conditions 
which may occur for the perturbed models, and in section |6| we list some relations between the 
Virasoro characters of related minimal models. Finally we give our conclusions in section |^. 

We give some details of the structure constants in the conformal field theories we use in 
appendix |^. 



2 Minimal models and their perturbations 



A conformal field theory is characterised by the central charge c of the Virasoro algebra, the 
sets of conformal weights {h,h} of the primary fields $a(z,z) ~ (j)h{z)(pji{z) and the three 
point couplings of these fields ($a|$f,(l)|$c)- The minimal models of conformal field theory 



were first described in |10|, and we denote them by Mr,s where r,s are coprime integers 



greater than 1. The Virasoro central charge in these models take the values 

,= l_fc££, (2,1) 

rs 

and in Air,s there are (r — l)(s — l)/2 possible values of h and h, given by 

(ms — nr)'^ — (r — s)^ 
h^n = i ^ , 0<m<r, 0<n<s, 2.2) 

although which actual pairs {h, h} are the conformal weights of fields is not determined a 
priori. For |r — s| = 1 Mr,s is unitary, and the lowest value of h is hi^i = 0; in all other cases 
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there is some allowed value of /i < 0, but there is a unified formula for the effective central 
charge c- 24hmin, 

6 

CeS = C - 2Ahmin = 1 ■ (2.3) 

rs 

One can consider perturbations of these conformally invariant field theories by the addition 
to the action of some term 



X j d'x (x) , (2.4) 

where ^(m,n) is the conformal primary field corresponding to h = h = hm,n- First-order 
perturbation theory calculations, backed up by arguments based on counting of states, suggest 
that the perturbed field theory is integrable for (m, n) one of (1, 2), (1, 3), (1, 5), (2, 1), (3, 1) 
and (5,1). For the unitary minimal models A^^.r+i, only (1,2), (1,3) and (2,1) are relevant 
perturbations, and consequently these have attracted the most attention to date. 

Much insight can be gained into minimal models and perturbed minimal models by con- 
sidering particular realisations as Lagrangian field theories. The Liouville action for a scalar 
field (j), 

So= f ( ^d,ct)d,(t> + exp(/3V20) \ , (2.5) 



. 47r 

provides a realisation of the minimal models with c = 13 -|- 6/3^ -|- 6/3~^. The minimal models 
correspond to the purely imaginary values of /3 = iyjr/s. Although this means that care must 



be taken with the interpretation of the action, Felder has provided a construction |11| which 



shows that the correlation functions defined in this theory are indeed those of the minimal 



models as found by Dotsenko and Fateev |12]: for the vacuum sector one must not consider 
the full Hilbert space of a scalar field but instead the reduced space Ker Q/Im : Q''^"^ : where 
Q is the screening charge 

Q = f dz exp(^/2</.//3) , (2.6) 



and : ^ : is an appropriately defined normal-ordered multiple integral. The vertex opera- 
tors 

exp{-am,n(p/V2) , am,n = (n- l)(3 - {m- l)/(3 , (2.7) 

transform as primary fields with conformal weight /im,n. Consequently, one can construct a 
formal perturbed conformal field theory action as 



S = So + \Jd^z exp{-am,n(l>/V2) . (2. 



For {m,n) = (1,3) this gives the standard sine-Gordon action, whereas for (m,n) = (1,2) or 
(1,5) we obtain the Bullough-Dodd or Zhiber-Mikhailov-Shabat (ZMS) model, which is the 

(2) 

affine Toda theory related to 02 . 

It is this latter model in which we are most interested, and to fix on a normalisation of 
the coupling constant, we shall say that the ZMS model with coupling constant 7 has the 
action 

SzMsil) 



J d^x (^{df^iff + ^ (exp(i v^(/p) + 2 exp(-i v^v?)) j , (2 
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where m is a couphng constant with dimension of mass. The connection between the variables 
in (|2.5D and (|2.9|) is given by 

4) = V^^ip , (3 = . (2.10) 

It is immediately apparent that the ZMS action may be viewed as a perturbed conformal 
field theory in two different ways: each of the two exponentials may be thought of as the 
perturbation and the remainder as the Liouville action: 

The actual values of r', s' depend on the arithmetic properties of r and s, since we require 
r', s' coprime. Hence the three distinct cases of table I can arise. 





r mod 4 


r' 


s' 


I 





r/4 


s 


II 


2 


r/2 


2s 


III 


1,3 


r 


4s 



Table 1: Relation between (r, s) and {r',s') 

These three cases differ in the way the effective central charges of Mr,s and Mr',s' are 
related. There are interesting examples of each of these cases as follows: 

I Here Ceg = 1 — 6/rs, Ceg' = 1 — 24/(rs) so that Cgg' < Ceg. Since the effective central 
charge is a measure of the number of degrees of freedom in the UV theory, we expect 
that there will be fewer particles in Mr',s' + *^(i,5) than in Air,s + ^(i,2)- 

The only (1,5) perturbation that has been investigated before falls into this class, being 
■M.2,9 + '^(i 5) which is related to A^g g + 2) and which was discussed by Martins et 
alin0,|.' 

II Here Cqs' = c^g and we expect the greatest similarity between the two models. 

We will discuss two models, A^6,5 + '^'(1,2) which is related to A^3,io + '^'(1,5)5 and + 
^(1,2) which is related to A^3,i4 + *^(i,5)- 

The first of these two models is unusual in that it can be expressed as a tensor product, 
M.3,10 = A^2,5 -A^2,5- The fields in M.3A0 can be classified into Z2-even and Z2- 
odd fields with respect to the Z2 map provided by flipping the tensor product. The 
perturbing operator $(1,5) is in the Z2-even sector and is nothing other than (li8)'5(i,2) + 
^(1,2) ® l)/\/2. As a result, 

-^3,10 + ^(1,5) = (-^2,5 + ^(1,2))^' , (2.12) 

is identically true in the Z2 even sector, and the spectrum consists of two self-conjugate 
particles of equal mass, each with the S'-matrix of the A^2,5 + '^(1,2) model. 
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In the second case the spectrum consists of six particles {/, L, h, h, H} for A^6,7 + ^(i,2) 
and {^4, A, C, B, B, D} for Ms^u + '^'(1,5) The S'-matrix in the second case is unusual in 
that it is not possible to define Sab and satisfying the usual unitarity condition, 
i.e. they are not pure phases. However, a formal application of the S-matrix-bootstrap 
gives Sab = Si^, = S'l'"^ as in (1:79), and the remaining 5-matrices can be found 
by application of the ^-matrix-bootstrap to the poles in Sab, S'^b- 

III In this case Ces' > Ces and we expect there to be more particles in M.r',s' + ^(1,5) than 

in Mr,s + ^{1,2)- 

Here we consider the case of the thermal perturbation of the Ising model, i.e. Ms^a + 
^(1,2) for which the corresponding model is Ms^iq + *^(i,5)- 

Although the bootstrap was not completed in [1], again we expect that A^3,i6 + *^(i,5) 
has the same peculiarities as Ms^u + ^(1,5) • 



3 TBA predictions from the 5*— matrices. 



As is well-known, the method of Thermodynamic Bethe Ansatz (TBA) is a powerful method 
for extracting the finite-size behaviour of a massive theory on a circle. One can also calculate 

this behaviour for the perturbed conformal field theory, both in perturbation theory about 
the conformal field theory and also numerically, and together these two methods provide a 
stringent test of the conjectured S'-matrices. 

We first formulate the conformal field theory on a cylinder of size R, and consider the 
Hamiltonian which propagates along the cylinder: 



H 

where 



HcFT + xjdx^z,z), (3.1) 



27r / - c 
HcFT = -f^iLo + Lo- — ] , (3.2) 



V 12, 

is the Hamiltonian of the ultraviolet CFT, ^{z,z) is some relevant operator corresponding 
to the perturbation and A is a dimensionful coupling. If the conformal dimension of the 
perturbing field $(z, z) is d$ = 2 — y, then A has mass dimension y. In terms of a specific 
mass unit m (which will be taken to be the mass of the lowest lying excitation of the model) 
we have 

A = amy , (3.3) 

where a is a dimensionless constant. The perturbing Hamiltonian is given by the second 
term: 

V = xjdx $(z, z) , (3.4) 
Perturbation theory in A gives the following expansion for the ground state energy 

n=l 
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The terms in the series come from the perturbative corrections, with coefficients 

(_-\\n j. 

Cn = ^—J-R'-''' / im n ^fe) )o.c , (3.6) 

where the integration is taken over the cyhnder and the indices 0, c refer to the connected 
contribution to the Green function in the original conformal field theory on the cylinder. For 
unitary theories, Ci is zero, while in the case of non-unitary models the following formula is 
obtained in [0], 

Ci = -(27r)-'^^^C<i,o$$o ' (3-7) 

where C^^^^g is the operator product coefficient for the perturbing operator $ and the 
operator with the lowest (negative) conformal weight $o- 

From the TBA, the same expansion can be computed using the conjectured scattering 
amplitudes. If we introduce the scaling length r = mR, where m is taken to be the mass of 
the lightest particle in the model, then the TBA yields Eq(R) as the series 



oo 



m 6r r 

n=l 



The first two terms, c and A, can be calculated analytically from the TBA whereas the a„ 
have to be calculated numerically. The agreement of the two series expansions (|3.5| ) and ( |3.8D 
leads to the following relations: 

Ceff = c, a„ = C„(mf A)" . (3.9) 



The linear term in (p.^), which is absent in ( |3.5| ), is interpreted as the bulk energy term in 
finite volume. 

The second relation in ( p.9| ) yields a mass gap formula in the special case n = 1, 

A = . (3.10) 

Cl 

The first nontrivial check is provided by the reproduction of the correct effective central 
charge, while the bulk energy constant and the mass-gap can be measured using TCSA. The 
n> 1 equalities could be used to perform further consistency checks, but we shall not do this. 

Now we give the TBA results for the models TWs^io and M^m by relating them to known 
results. 



3.1 TBA for Mzin + ^ 



3,10 y^(l,5) 



It was pointed out in |18] that the TBA calculation for the scaling Yang-Lee model (which is 
-^2,5 + '^'(1,2)) is identical to the TBA calculation for the scaling Potts model M^fi + '^(2,1) 
for the following reason: The scaling Yang-Lee model contains one scalar particle B with the 
5-matrix 

Sbb={\){1). (3.11) 
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while the spectrum of the scaHng Potts model consists of a conjugate particle pair A and A 
with the S'-matrix 

Saa = (^) , = [I) ■ (3.12) 

The two S'-matrices satisfy 

SaaS^a = ^BB ■ (3.13) 

As shown in |18], this means that the TBA equation in the charge symmetric Gibbs state for 
the scaling Potts model is equivalent to two copies of the scaling Yang-Lee model. 

However, the model A^3,io + '^'(1,5) is explicitly equal to two copies of the Yang-Lee model, 
and therefore the ground state energy computed from the TBA of the model Ai^fi + ^(2,1) is 
identical to that of the model A^s^io + ^(1,5)- 



3.2 TBA for 7^3,14 + ^>(i,5) 

The model A^3,i4 + $(i,5) is related to A^6,7 + *^'{i,2)- The S'-matrix of the latter is the so-called 
minimal eg S'-matrix, derived in . If we denote the S'-matrix of the i and j particles in 



the non-unitary model by S^j'^^^ , while that of the corresponding particles in the unitary 
case by Sj-^'^\ then we have the following interesting relation, 

^.(3,14) ^.(3,14) _ ^.(6,7) ^{6,7) , . x 

even for the case when S'|^'^^ is not a pure phase. 

Sab Sy^j^ = Sih Sij^ . (3.15) 

This implies that the ground state energy calculated using the charge symmetric Gibbs state 
for TBA is the same in the two models. 



The TBA for + *^(i,2) has been analysed in |17], where they find 

c = 6/7 , A^^^ = , ai = 0.0027765 . . . , (3.16) 

6 -|- 2v3 

For A^3^i4, c = —114/7 and the operator $0 with the minimal dimension is ^*(i,5) with 
hmin = —5/7, so that Ceff = 6/7 as predicted. Using d^^^^ = —20.3634. . ., the mass-gap m 
is given by 

\ = an?^/\ a = -0.011833... . (3.17) 



3.3 TBA for 7^3,16 + ^(1,5) 

The arguments used in the two previous subsections to show that the TBA analysis of the 
two related models are equal are not valid here since the non-unitary model contains more 
particles in its spectrum than the corresponding unitary one. Therefore one does not expect 
the TBA prediction derived from the Ising model results to hold in this case. 
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For example, consider the effective central charge: In the Ai^^w + ^(1,5) model we have 
c = —161/8, hmin = —7/8 (corresponding to the operator <&(i,5)) and so Cefr = 7/8. The 
related model is ^^3,4 + ^(1,2) with Cefj = 1/2, and which has 8 particles with the minimal eg 
5-matrix, for which c = 1/2 as required. 



In the eg TBA [17], the mass gap is expressed in terms of the lowest lying breather. If we 
try to use the eg TBA expansion for A^3,i6 given in [^], we get 

ai = -0.001236... , C$$o<i, = 51.85 ... , (3.18) 

and therefore the predicted mass gap is 

A = am^^/^ , a = -.003735 . . . (3.19) 

The bulk energy constant can be obtained from |17] as follows: 

A^^"^ = ^ = .06173 . . . (3.20) 

4 [sin + sin f + sin ^ + sin ^ + sin ^ + sin ^ j 

We give these simply to show that they do not agree with the numerical results of the TCSA 
in section ^ as indeed they should not. 



3.4 BA predictions from the 5 matrix 

The TBA system is a set of complicated non-linear integral equations for the ground-state 
energy Eq{R), and even more complicated equations for higher energy levels There is a 
much simpler set of equations which give the dominant behaviour of two-particle excitations, 
which comes simply from the usual Bethe Ansatz for a two-particle wave function as follows. 

On a circle the momenta of two-particle states are quantised. If we consider the zero 
momentum sector with particles i and j of masses mj, mj, momenta p, —p, and total energy 
E{R), then we have 

p = mi sinh(6'i) = -ruj sinh(6'2) , E{R) - Eo{R) = \J p^ + mf + ,Jp'^ + mj , (3.21) 

where Eo{R) is the ground state energy. If the particles have purely elastic scattering with 
5-matrix Sij{9) then the quantisation condition is 

expiipR) Sy{di-e2) = ±1, (3.22) 

with a —1 for i and j both fermions and 1 otherwise. Consequently it is possible to extract 
the S-matrix, or more properly the phase shift 5{9) = —ilog{S{9)), from the function E{R) — 
Eq (R) . If there are two or more particles of the same mass with non-diagonal scattering then 
the quantisation condition is ( |3.22| ) with Sij{6) replaced by the eigenvalues of the two-particle 
S'-matrix. We are able to identify several such two-particle lines in our numerical results and 
we can compare our conjectured S-matrices with those of the related models. 

However, there is an ambiguity in our measurements of S{6) which we should mention. 
The particle statistics are determined by the sign ±1 in equation ( |3.22| ); there is also the 
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notion of 'type' which enters the TBA and which is given by (±1)5(0). The folklore is that 
all particles in theories such as ours are of fermionic type, that is (±1)S'(0) = —1, and this is 
the convention we have used in our TBA and BA calculations. However, the TBA and BA 
predictions are unaltered by S{6) — > —S{9) and fermion^boson, and this ambiguity can only 
be fixed by requiring the correct signs in the S-matrix bootstrap. 



3.5 Extension of BA predictions from the S matrix. 

There are two situations in which we would like to consider extending equations ( |3.21| ) and 
( 3.22| ) to complex rapidities. 



Firstly, if the S'-matrix is unitary, i.e. a pure phase for real rapidity, then we can reasonably 
expect E{R) to be real. This requires that either 9i are real and p is real, or 9i is purely 
imaginary and and p is purely imaginary. In the second case, the energy will be less than 
the sum of the particle masses, and so cannot correspond to an excited state of these two 
particles. However, if one considers the full TBA equations for the single particle state fc, 
then it is quite possible that for large R that they are dominated by this BA equation, if k 
appears as a bound state pole in Sij. Although this may give the leading corrections for large 
i?, as is the case for the Yang-Lee modelQ to the M.2,'0 + *^(i,2)) typically this will not be a 
good approximation for small R as other corrections become dominant. However, we shall 
see that this need not be the case. 

Secondly, if the ^-matrix is not a pure phase for real rapidity difference, then requiring 
E{R) real will lead to no solution of ( |3.22| ). However, in this case, we can instead demand that 
R is real, which is certainly a requirement, and which will typically lead to complex values 
of 9i and E{R). We have the first such occurrence in + *l'(i,5)) where Sab and 

are not pure phases. They are exchanged under crossing, and as a result, the values of E{R) 
which we find by applying the BA to these two S'-matrices are conjugate. The BA is derived 
under the assumption that the particles are well separated and with real momenta, which is 
certainly not the case here; however it is remarkable that this simple rule does indeed describe 



the spectrum very well, and we compare these values with the TCSA data in subsection 4.1 



4 Numerical analysis in the Truncated Conformal Space Approach 

The Truncated Conformal Space Approach (TCSA) consists of calculating the matrix ele- 
ments of ([4.ID exactly for a finite dimensional subspace 7i of the full Hilbert space, and then 
numerically diagonalising the resulting matrix, as explained in [^]. Following the discus- 
sion in [^], the spectrum of the truncated Hamiltonian hxcsA exhibits three distinct scaling 
regimes; r small in which the UV conformal field theory dominates, r large in which the 
truncation effects dominate and an intermediate region in which one hopes that the spectrum 
will approximate that of the massive theory on a circle. The TCSA method is fully described 

^One can check that the leading correction to the single-particle mass found in is exactly given by the 
BA equation at imaginary rapidity 
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in e.g. |^0|, but we recall the main elements here. The full Hamiltonian is: 



27r / ry \ 

hir) = -{H, + aj-^H,) , (4.1) 

in units in which m = 1 and where 

Ho = Lo + Lo-^. (4.2) 

The matrix elements of the perturbing Hamiltonian can be evaluated using the formula 

27r Ry 

{"fflVm = -J^^^^i'^fimmpUne , (4.3) 

where the matrix element of the perturbing field is taken on the conformal plane and evaluated 
using standard conformal field theory techniques. From this the operator Hi can be identified 
with the perturbing field ^(1) itself. 

There are some symmetries which make the calculation of hxcsA simpler: firstly, since the 
Hamiltonian is Lorentz invariant, we can restrict attention to the spin-zero sector of the full 
Hilbert space. Secondly, in our models the perturbing operator $(1) = $(1,5) (1) preserves 
distinct sectors of the full Hilbert space, and as a result we may diagonalise hxcsA on these 
sectors separately. We can also choose to normalise the highest weight states as 

= ±5,„ (4.4) 

where the signs are chosen to ensure the reality of the three-point couplings, which are them- 
selves given in appendix 

As discussed in section ^, the following numerical results can be extracted from the TCSA: 
(1) The ground-state energy: in the scaling regime the ground state energy should be ap- 
proximately £'o(r)/m ~ Ar, and one can measure A. (2) The mass-gap can be found and 
compared with the TBA predictions. In practice, we adopt the value of a from the TBA and 
then we expect that the lightest particle has mass 1. (3) Some of the two-particle lines can 
be identified, and for those, the two-particle energies can be compared with the Bethe Ansatz 
results derived from the ^-matrix. We now present the results of the Truncated Conformal 
Space calculations for the two models ^^3,14 + ^{i^b) ™d A^s^ig + *^(i,5)- 



4.1 The TCSA results for A43,i4 + $(15) 

The allowed values of h are hi^i, < ? < 14. It is possible to consider several different 
sectors in which ^*(i,5) acts locally but which are not themselves necessarily mutually local. 
In particular, we can consider the three sectors 'even', 'odd' and 'twisted' for which the fields 
all have integer spin, with the following sets of /i-ij} = [i, j], 

Sector Field content 

even [1, 1], [3, 3], [5, 5], [7, 7], [9, 9], [11, 11], [13, 13] 

odd [2,2], [4,4], [6,6], [8,8], [10,10], [12,12] ^ ' > 

twisted [1,13], [3,11], [5,9], [7,7], [9,5], [11,3], [13,1] 
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The coupling constants for first two sectors are those of the 'A' type modular invariant, which 
we give in appendix |A|, the even and twisted sectors correspond to the 'D' type modular 
invariant. 

To establish that this model has the spectrum we have predicted, it is only necessary to 
consider the 'even' sector. We truncated the space by requiring Lq + Lq < 10 which leaves 
285 states. In figure || we give the real eigenvalues of h{r) and pick out the first 13 lines; it is 
clear that scaling has set in for the first three eigenvalues by r = 10. 



E [r] r E [r] 




(a) Plot of E{r) versus r. (b) Plot of rE[r) versus r. 



Figure 1: The first 13 eigenvalues of h{r) in A^3,i4 + ^(1,5) • 

Some of the lines appear to terminate - in fact two real eigenvalues collide and become a 
complex conjugate pair of eigenvalues; this is due to truncation effects as explained in |^]. 

From figure ||a it is possible to estimate the slope of the ground-state energy and the first 
three excitation energies in the scaling region (which we take to be 13 < r < 17), and we give 
these in table § for truncation levels 8, 9y and 10 to show systematic errors together with 
estimates of the error of the measurements, together with the values predicted by the TBA 
and perturbation theory. 



Level 


dim'H 


A 


mi 


m2 


m3 


8 


136 


-0.105 ±0.05 


0.994 ± 0.006 


1.420 ±0.006 


1.94 ±0.03 


gi 
^7 


280 


-0.106 ±0.02 


1.01 ±0.02 


1.417 ±0.002 


1.95 ±0.03 


10 


285 


-0.105 ±0.05 


0.996 ± 0.004 


1.415 ±0.004 


1.95 ±0.03 


Exact 




-0.1056... 


1 


1.414... 


1.932... 



Table 2: Results for A^s^u ± ^{1,^) 
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However, as we discuss further in section ^, the particles A and A, and B, B we used to 
diagonaUse the S-matrix in section Q are not invariant under Z2 and consequently there is 
no direct connection between these particles and the excitations in the even sector; rather the 
particle of lowest mass in the even sector is a symmetric combination of A and A; similarly, 
the zero-momentum two-particle states of A and A in the even sector are \K(6)K{—6)) and 
\K{6)K{—9)). However, it is still true that the eigenvalues of the S-matrix on these two 
particle states are the matrix elements Saa{G) and 5'^^(^), and we have indicated the first 
two lines corresponding to the eigenvalues by dotted and dashed lines respectively. 

From figure [l|b we can extract the S'-matrix element Saa{(^), which we compare with the 
conjectured results in figure the data is from the second energy level which corresponds 
to these particles, i.e. the 10th eigenvalues of h{r) for small rQ. As is usual, there is a small 
discrepancy for small r between the Bethe-Ansatz predictions for E{r) and the actual values 
due to finite-size effects, but otherwise the agreement is very good. In figure ^3 we give the 
same data for 5cc(^) which corresponds to the 13th eigenvalue of h{r) for small r. The 
results for Sbb and S"^^ are equally convincing. 



I Log[S] 

7r 




(a) The AA S'-matrices. (b) The CC 5-matrix. 



Figure 2: The S'-matrices as extracted from the TCSA compared with theory in A^3^i4+$(i 5). 



The interpretation of the lowest energy solution of equation (|3.22 ) for S = (the 5th 



eigenvalue for small r) is extremely interesting. Since Saa{0) = — 1, the particles A and A 
are bosons, and so the Bethe Ansatz for the nth AA excitation is 

r sinh(6l) + 6^^(29) = 2mr . (4.6) 

Since Saa{^) —0{3 + 2\/2) for small 9 we see that equation ( [4. 61) with n = only has 
solutions for real 9 for r < 6 + 4\/3 = 12.92 .... However, for r large, we have identified this 

2 In figs. |[i and ||b, small r corresponds to large 6, and large r to small 6. 
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5th line as the single particle Z2 even component of the higher kink doublet B, B, and hence 
obtained the estimate of the mass = 1.95 . . .. 

We can reconcile these two interpretations when we notice that for r > 6 + 4\/3 the solution 
to equation ( [4. 61) corresponds to 9 imaginary. As r — > 00, we see that 9 ivr/G, and there is 
a pole in S^^{9) at 9 = in/G which corresponds to the particles B and B. Thus we see direct 
evidence that the first excited kink (corresponding to the particles B and B) appears as the 
direct channel pole of KK at 9 = in/ 6. 

As a consequence, the TCSA gives the value of the 5-matrix Sj^^{9) for all real 9 and also 
for the imaginary values between and z7r/6. In figure ^ we plot the theoretical and numerical 
values of both —6j^^{9) versus 9 and log{S{9)) versus —i9. There is excellent agreement for 
all real values of 9, and also for imaginary values of up to r ~ 17 (corresponding to ~ OAi) 
where scaling starts to break down and truncation effects take over^ 



I Log[S] Log[S] 




(a) Sj^^ for real 9. (b) Sj^^ for imaginary 9. 

Figure 3: The S'-matrices as extracted from the TCSA compared with theory in A^s, 14+^(1, 5)- 

However, figure || does not tell the whole story, as we have only plotted the real eigenvalues. 
There are also eigenvalues which are genuinely complex, and in figure ^ we show the real and 
imaginary parts of those eigenvalues A amongst the first 60, for which 9'm(A) > 10~^. While 
some of these are simply the results of truncation error, there are also several series of lines 
which are not, and for these we have included the predictions of the 'two-particle BA' using 
the ^-matrices Sab and and as can be seen, there is excellent agreement. 



■^In figs, bk and Hb, small r corresponds to large 9, and r— ^ootoS— s-ivr/G 
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Re [E [R] ] 



(a) 5Re(A). 



Im[E [R] ] 




0.02 



-0.02 



-0 . 04 




2.5 5 7.5 10 12.5 15 17.5 20 



,.„, .•..«ipl";"»t 



(b) 9m(A). 



Figure 4: The complex eigenvalues h{r) in 7^3,14 + '^'(1,5): data and BA predictions 



4.2 The TCSA results for Ms^m + $(1,5) 



The allowed values of h are hi^i, < i < 16. It is again possible to consider several different 
sectors in which ^*(i,5) acts locally but which are not themselves necessarily mutually local. 
In particular, we can again define three sectors 'even', 'odd' and 'twisted' for which the fields 
all have integer spin, with the following sets of {hi^i,hij} = 



Sector Field content 

even [1, 1], [3, 3], [5, 5], [7, 7], [9, 9], [11, 11], [13, 13], [15, 15] 

odd [2,2], [4,4], [6,6], [8,8], [10,10], [12,12], [14,14] ^^'^^ 

twisted [2, 14] , [4, 12] , [6, 10] , [8, 8] , [8, 8] , [10, 6] , [12, 4] , [14, 2] . 

The coupling constants for first two sectors are those of the 'A' type modular invariant, which 
we give in appendix the twisted sector corresponds to the 'D' type modular invariant. 

To establish that this model has the spectrum we have predicted, it is again only necessary 
to consider the 'even' sector. We truncated the space by requiring Lq + Lq < 12 which leaves 
870 states. In figure |5| we give the real eigenvalues of h{r) with the first 14 lines picked out. 
It is less clear that scaling has set in for the first five eigenvalues, and both the masses of 
the single particles and the S'-matrices as extracted from the excited two-particle lines are 
in good agreement with theory. However, to achieve a reliable scaling regime it has been 
necessary to use a much larger truncated space for this model than the previous one; indeed 
when the space is truncated to 285 states the lines corresponding to the third particle are no 
longer present, and for that reason we leave the measurement of its mass blank in table |^. 
The severity of the truncation effects which still remain can also be seen in the large gaps 
in the 2nd, 3rd, 4th and 5th lines (when the eigenvalues become complex) and which would 
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presumably disappear when the truncation level is increased. 

E [r] r E [r] 




(a) Plot of E(r) versus r. (b) Plot of rE[r) versus r. 

Figure 5: The first 14 eigenvalues of h{r) in A^3,i6 + *^'(i,5)- 

From figure ||a it is possible to estimate the slope of the ground-state energy and the first 
three excitation energies in the scaling region (which we take to be 10 < r < 12), and we give 
these in table ^ for the truncation levels 8|, 9| and 12 to show systematic errors together 
with estimates of the error of the measurements. 



Level 


dim'H 


A 


mi 


m2/mi 


m^/mi 


8| 


285 


-0.08 ± 0.02 


1.01 ±0.06 


1.14 ±0.04 




9! 


369 


-0.081 ±0.05 


0.96 ±0.03 


1.175 ±0.02 


1.81 ±0.02 


12 


870 


-0.0815 ±0.0008 


0.975 ± 0.006 


1.178 ±0.02 


1.83 ±0.02 


Exact 








1.175... 


1.827... 



Table 3: Results for Mz,i6 + 

It is clear from the measurement of A that this is not given by the standard minimal 
eg TBA, as should be the case, and that the mass ratios are consistent with those of the 
S'-matrices of |||]. 

From figure |5|b we can extract the eigenvalues of the S-matrices of the fundamental kinks, 
which we compare with the conjectured results in figure ^. The data points come from the 
9th, 11th and 23rd eigenvalues of h{r) (for r small) with the continuous lines being Si{9), 
S2{0), and the dashed lines S3{9) and S4{9). 

We can also identify the first two-breather line, which is the 14th eigenvalue of h{r) for small 
R, and we compare the S-matrix extracted from these numerical results with the theoretical 



15 



result of 



2m sin ■ 



vr 
15 



1.175. . .m ,SbiBi 



(4.8) 



in figure |6|. In this case we have SsiBii^) = 1 so that we take Bi fermionic for the Bethe 
Ansatz, in accord with the folklore that bosonic and fermionic particles have ^(O) = — 1 and 
1 respectively; in any case there is agreement between theory and experiment for this choice. 



I Log[S] 




I Log[S] 

lOr 



. Theta 



(a) The Kink-Kink S-matrix eigenvalues. (b) The Bi-Bi S-matrix. 



Figure 6: The 5-matrices as extracted from the TCSA compared with theory in A^3,i6 + *^(i,5) 



5 Z2 symmetry and sectors 



Minimal models of CFT are known to exhibit Z2 symmetry, which has been used above to 
classify the sectors of the Hilbert-space. These sectors are preserved under $(1,5) perturbation 
and therefore it is natural to expect a corresponding symmetry of the massive field theory. 

There is a natural Z2 action on the particles in the unrestricted ZMS model, completely 
analogous to the classification of the irreducible representations of the classical group SU{2) 
by its centre Z2 , with integer spin representations corresponding to the trivial (even) class and 
half-integer spin representations to the nontrivial (odd) class. Although the Uq{sl{2)) action 
does not survive RSOS restriction, the Z2 action does, and as a consequence, the kinks K are 
even and K odd. Since the Z2 symmetry is a remnant of the full quantum group symmetry, 
it commutes with the 5-matrix, and it is possible to simultaneously diagonalise both the 
S'~matrix and Z2 on the two-particle states. However, it is not possible to do this on one- 
particle states. As a result, the one-particle states which diagonalise the 5-matrix (denoted 
A, A in case of A^3,i4 and A^s^ie and A, B and the case of ^^3,10) are not Z2 eigenstates; 
rather, they are exchanged with each other under Z2. Therefore the natural description of 
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the content of each sector is m terms of the states K and K, which are self-conjugate and 
Z2-even or odd, respectively. 

We now examine the content of the three different sectors (even, odd and twisted) of the 
simplest model A^3,io + *^'(i,5)i as an example. Exactly as for the models Ais^u and M-s^iq 
the values of = of A^s^io fall into even, odd and twisted as follows: 

Sector Field content 

even [1, 1], [3, 3], [5, 5], [7, 7], [9, 9], 

odd [2,2], [4,4], [6,6], [8,8], ^^'"^ 

twisted [1,9], [3,7], [5,5], [7,3], [9,1]. 

The direct sum of the odd and twisted sectors correspond to the product (A^2,5 + ^(1,2))*^^) 
and give the D type modular invariant in the conformal limit. The direct sum of the even and 
odd sectors gives the A-type modular invariant partition function. Direct TCSA investigation 
of the odd sector gives very similar pictures to those of the even sector in the scaling regime: 
both sectors have a ground state, a one-particle line and a very similar structure of two- 
particle lines. The degeneracy of the ground state energies is broken by an amount which is 
exponentially small for large r. The picture is very similar for ^^3.14 -|- ^(1,5)- 

For the model A^(3,i6)) the four eigenvalues of the kink-kink S'-matrix are different, so 
that we now expect that the two-particle lines in the even and odd sectors are given by 
different Bethe Ansatz results. In the even sector, we expect the two-particle lines to be 
linear combinations of the states \K{6)K(—0)) and \K{6)K(—6)), whereas those in the even 
sector to be linear combinations of \K{9)K{—0)) and \K{9)K(—9)) . The identification above 
is consistent with the fact that the eigenvalues on the Z2 even sector, {S3 and S4, the dashed 
lines) give a better approximation to the phase-shift extracted from the even sector, than 
those on the Z2 odd sector {Si and 5*2, the solid lines). 

This picture can be formulated generically as follows: states corresponding to RSOS se- 
quences {jm ■ ■ ■ j2, ji} which start and end with an even vacuum state (i.e. ji and jn are 
integer-spin representations, which in our three cases means the state |0)) originate from the 
even sector of the underlying CFT, and the ones starting with an odd vacuum and ending 
with an odd vacuum derive from the odd sector. 



6 Character Identities 



As has been shown extensively, there are interesting sum forms for Virasoro characters which 



may be related to the particle structure of perturbed minimal models (see e.g. [21] for the 
cases related to the minimal eg 5-matrix.). Since the breather sectors of our related models 
are the same, and the characters of the Virasoro representations in these models may be given 
in terms of sum formulae, we might expect to find relations between the Virasoro characters 
in the two models. This is indeed the case, as we show now. The character of the irreducible 
representation Lc^h of the Virasoro algebra of central charge c with highest weight h is defined 
as 
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For pairs Air,s and J^r',s' of type II (r' = r/2, s' = 2s) we find that 

X.Cr,s,h^/2 — 2m.n Xcr.s ,h,. /2 — 2m,s — n -^'^r' , s' ''*m,3'/2- 2n -^"^r' ,s' I'^r' -m.s' /2 — 2n ' 

Xcr,s,h2?n.n ~ Xc^l ,.1 ,hm,2n' (6-1) 

Note tliat these identities do not involve all the characters of either model. For the pair A^s^io 
and M.^fi these relations give 

Xc3,io,/ll,l ~ Xc3,io,/ll,9 Xc5, 6,^2,1 ~ Xctj^Q,hz,l 1 

'^C3,10,/ll,3 ~ Xc3_io,fel,7 Xcs^Bi'^l,! ~ Xc5_6,h4_i ) 

Xc^,lQ,h-i^2 Xc5,6,fci,2 ; 

Xc3,l0,fel,4 ~ Xc5,6,/i2,2 ' 

/^C3, 10,^1, 6 Xc5,6,/l3,2 ; 

Xc3,l0,/ll,8 ~ Xc5,6,ft.4,2 • (6-2) 

For the pairs A^r,s and Mr,As of type III we find that 

_ f Xcr,4s,ft.{r-m)/2,2{a-n) ~ '*(ICr,4s ,/l(r-m)/2,2(s + n) "-"^d 
XCr s^hm,n I V i. — V L TTI PVPTl * \ * / 

I, ACr,4s,Al„/2,2n AC,.,4s Am/2,4s- 2n //tCVCil 

In the case of Ms^iq and 7^3,4 this yields 

Xc3,4,/ii,i Xc3,i6,/ll,6 ~ Xc3_i6,/ll,l0 ) 

Xc3,4,/il,2 ~ Xc3,i6,/ll,4 ~ Xc3, 16,^1, 12 ' 
Xc3,4,/ll,3 — Xc3,i6,/ll,2 ~ Xc3,i6,/ll,14 ' 

(6.4) 



Analogous relations for type I pairs are obtained by swapping r s',s ^ r' and m 



^ n. 



All these relations may be proven easily using the formulae in [22|. Whether the identities 
in equations ( |6.lD , (|6.3| ) have any real interpretation in terms of the particle structures of 
PCFTs is an interesting question to which we hope to have an answer shortly, although it is 
interesting to note that in each case we have investigated when the difference of two characters 
appears, it still has an expansion in q with positive coefficients. 



7 Conclusions 



The above results make it possible to clarify the connection between $(1,2) and '^(1,5) per- 
turbations of minimal models. We have shown that they can be related through the ZMS 
Lagrangian and that the relation falls into one of the three possible classes of type I, II or 
III, according to the relative magnitude of the effective central charges. Using the S-matrices 
obtained via RSOS restriction in [Q], we tested them against the predictions of the TCSA 
and TBA method. We have confirmed the picture of *I'(i,5) perturbations based on the RSOS 
restriction and have shown that the two related theories are different reductions of the same 
underlying affine Toda field theory. The existence of such different restrictions relies upon the 
existence of non-equivalent sl{2) subalgebras of the affine Kac-Moody algebra 02 underlying 
the ZMS model. It appears reasonable to expect that such relations should exist between 
perturbations of Pi^-minimal models corresponding to other imaginary coupling affine Toda 
field theories. 
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The results show that in the type II case (models A4(3 iq) and A^(3,i4)) we end up with 
the same mass spectra, although different scattering theories, whereas in the cases of type I 
(of which the example of M.2,9 was treated in ^) and type III (of which A^s^ie is an 
example) even the mass spectra are different just as expected. In the two type II cases the 
TBA equations turned out to be the same as those of the corresponding <I>(^ 2) "Perturbed 
conformal field theories. 

We have also shown that the strange features of the S-matrices for 7^3,14 + ^(1,5) found 
in Q are in fact correct. This model does indeed break unitarity much more severely than 
the '^(1,2) ^-iid $(1,3) perturbations considered before. The spectrum is not entirely real, and 
it is well described by the S'-matrices no longer being pure phases. We believe that this is a 
generic feature of RSOS restrictions, and that the previous situation with a real spectrum is 
in fact rather exceptional. Indeed, looking at the 5-matrices found for Toda theory in 



1 24 1 we see that the S-matrices for soliton-breather scattering (which should not be altered 



by RSOS reduction) are also not pure phases in general. 

We have also shown that in one case, the Bethe Ansatz equations give a good approximation 
for a single-particle state energy by continuation to imaginary rapidity. While it is easy to 
check that in general this is not the case, it is true that the leading corrections to the particle 
mass in A^2,5 + ^(1,2) filso given by the BA equations, and that this can be seen by 
examining the large R behaviour of the full TBA equations describing this energy Q. We 
believe that in our case, the full TBA equations are simply dominated by the BA equations 
for all values of R, rather than simply for large R. 

Unfortunately, due to our lack of knowledge about the conformal field theory three-point 
coupling constants for D-type modular invariants of non-unitary models, we have not been 
able to extract precise TCSA predictions for the twisted sector, although preliminary results 
(using some guesses for the structure constants) have turned out to be in good agreement 
with our expectations based on the Z2 picture. It seems worthwhile to tackle the problem of 
the coupling constants in the non-unitary minimal models and in this way to facilitate the 
further investigation of their structure as well. 

Since the RSOS restriction does not modify the breather sector of the affine Toda field 
theory, one can expect that this sector is identical in corresponding '&(i,2) and '^(1,5) pertur- 
bations irrespective of the type I, II or III nature of the pair. This leads to the existence of 
'strange' character identities between characters of different minimal models. In addition to 
that, in the type II case we expect relations between the states corresponding to the solitonic 
sectors as well. Some character identities of this type have been described, although it is not 
clear whether they originate from some connections between different PCFTs in the case of 
general minimal model pairs. 
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A Coupling constants in non-unitary conformal field theories 



The coupling constants for fields of type ^{i,a) &re given in |12| in an asymmetric fashion: 
they have two sets of fields fa, 'fa which for our purposes are just (different) multiples of 
^(i,a) Using the results of |12] we arrive at 

{fa\ fbi'i^) \ fc) = Cabc , {fa\ fb) = C'aW 

where the coupling constants are 



b+c—a—l • 



a+c—b—l ' 



a+b—c—l ■ 



r]{ 



a+b+c-l ' 



abc 



e(a-l) ^b-l) e(c-l) 



rj{l) 



where 



(A.l) 



C(a) = n ^(/'•?') ' ^(") = n ~ ^) ' ^(^) = r(x)/r(l - x) , p = r/s . 

3=1 j=i 

From the unnormalised coupling constants Cabc it is then easy to choose normalised fields 



(l,a) 



^ifa such that ($(l,a) 



<f(l,6)(l) 



^(i,c)) is real and ($(i,a) 



(l,a) 



±1. 



We do not give the coupling constants for the D type modular invariant involving the fields 
in the twisted sector as these are not to be found in |^]. Petkova and Zuber ||2^ give these for 
the unitary minimal models, and it would be interesting to extend their results to all minimal 
models, but we do not do that here. 



References 



[1] G. Takacs, A new RSOS restriction of the Zhiber-Mikhailov-Shabat Model and ^(1,5) 



perturbations of nonunitary minimal model, preprint DAMTP-96-27, tLep-th/9604098 



[2] T. HoUowood and P. Mansfield, Rational Conformal Field Theory at, and away from criticality, 
as Toda Field Theories, Phys. Lett. B226 (1989) 73-. 

[3] D. Bernard and A. LeClair, Quantum Group Symmetries and non-local Currents in 2D QFT, 
Commun. Math. Phys. 142 (1991) 99-138. 

[4] F. A. Smirnov, Exact S -matrices for (f)i^2 perturbed minimal models of conformal field theory, 
Int. J. Mod. Phys. A6 No. 8. (1991) 1407-1428. 

[5] A. Koubek, S -matrices of (j)i^2~perturbed minimal models, IRF formulation and bootstrap 
program. Int. J. Mod. Phys. A9 No. 11 (1994) 1909-1927, |hep-th/9211134 



20 



A. Koubek, M. J. Martins and G. Mussardo, Scattering matrices for 0i,2 perturbed conformal 
minimal models in absence of kink states, Nucl. Phys. B368 (1992) 591-610. 

M. J. Martins, Constructing an S -matrix from the truncated conformal approach data, Phys. 
Lett. B262 (1991) 39-42. 

P.E. Dorey and R. Tateo, Excited states by analytic continuation of TEA equations, Durham 



preprint DTP-96/29, hep-th/9607167 



V.V. Bazhanov, S.L. Lukyanov and A.B. Zamolodchikov, Integrahle quantum field theories in 



finite volume: excited state energies. Preprint CLNS 96/1416, hep-th/9607099 



V. P. Yurov and AL B. Zamolodchikov, Truncated conformal space approach to scaling 
Lee-Yang model. Int. J. Mod. Phys. A5 No. 16 (1990) 3221-3245. 

A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, Infinite conformal symmetry in 
two-dimensional quantum field theory, Nucl. Phys. B241 (1984) 333-380. 

G. Felder, ERST Approach to Minimal Models, Nucl. Phys. B317 (1989) 215-236. 

VI. S. Dotsenko and V. A. Fateev, Conformal algebra and multipoint correlation functions in 
2D statistical models, Nucl. Phys. B240 [FS 12] (1984) 312-348; Four-point correlation 
functions and the operator product algebra in 2D conformal invariant theories with central 
charge C <l, Nucl. Phys. B251 [FS 13] (1985) 691-734. 

T. Eguchi and S.-K. Yang, Sine-Gordon theory at rational values of the coupling constant and 
minimal conformal models, Phys. Lett. B235 (1990) 282-286. 

G. Felder and A. LeClair, Restricted quantum affine symmetry of perturbed minimal conformal 



models. Int. J. Mod. Phys. A7 (Proc. Suppl.) 239-278, |hep-th/9109019| ; 

V. Pasquier and H. Saleur, Common structures between finite systems and conformal field 
theories through quantum groups, Nucl. Phys. B330 (1990) 523-556. 

F. A. Smirnov, The perturbated c < 1 conformal field theories as reductions of sine- Gordon 
models. Int. J. Mod. Phys. A4 No. 16 (1989) 4213-4220. 

A. LeClair, Restricted sine-Cordon theory and the minimal conformal series, Phys. Lett. B320 
(1989) 103-107. 

T. R. Klassen and E. Melzer, The thermodynamics of purely elastic scattering theories and 
conformal perturbation theory, Nucl. Phys. B350 (1991) 635-689. 

Al. B. Zamolodchikov, Thermodynamic Eethe ansatz in relativistic models, scaling 3-state Potts 
and Lee-Yang models, Nucl. Phys. B342 (1990) 695-720. 

G. Sotkov and C. J. Zhu, Eootstrap fusions and tricritical Potts model away from criticality, 
Phys. Lett. B229 (1989) 391-397. 

M. Lassig, G. Mussardo and J. L. Cardy, The scaling region of the tricritical Ising model in two 
dimensions, Nucl. Phys. B348 (1991) 591-618. 

A. Koubek, A method to determine the operator content of perturbed conformal field theories, 



Phys. Lett. B346 (1995) 275-283, tiep-th/9501028 



M. Rocha-Caridi, Vacuum vector representations of the Virasoro algebra, in 'Vertex Operators 
in Mathematics and Physics', MRSl publication # 3, Springer (1984), 451-473. 



21 



[23] V. Petkova and J.-B. Zubcr, On structure constants of sl{2) theories, Nucl. Phys. B438 (1995) 
347-372, |hep-th/9410209| . 

[24] G. M. Gandenberger, Exact S-matrices for bound states of a^"^ affine Toda solitons, 
Nucl. Phys. 449 (1995) 375-405, ^lep-th/9501136|. 



22 



